Effects of selection for ethanol preference on gene expression in the nucleus accumbens of HS-CC mice.

Previous studies on changes in murine brain gene expression associated with the selection for ethanol preference have used F2 intercross or heterogeneous stock (HS) founders, derived from standard laboratory strains. However, these populations represent only a small proportion of the genetic variance available in Mus musculus. To investigate a wider range of genetic diversity, we selected mice for ethanol preference using an HS derived from the eight strains of the collaborative cross. These HS mice were selectively bred (four generations) for high and low ethanol preference. The nucleus accumbens shell of naive S4 mice was interrogated using RNA sequencing (RNA-Seq). Gene networks were constructed using the weighted gene coexpression network analysis assessing both coexpression and cosplicing. Selection targeted one of the network coexpression modules (greenyellow) that was significantly enriched in genes associated with receptor signaling activity including Chrna7, Grin2a, Htr2a and Oprd1. Connectivity in the module as measured by changes in the hub nodes was significantly reduced in the low preference line. Of particular interest was the observation that selection had marked effects on a large number of cell adhesion molecules, including cadherins and protocadherins. In addition, the coexpression data showed that selection had marked effects on long non-coding RNA hub nodes. Analysis of the cosplicing network data showed a significant effect of selection on a large cluster of Ras GTPase-binding genes including Cdkl5, Cyfip1, Ndrg1, Sod1 and Stxbp5. These data in part support the earlier observation that preference is linked to Ras/Mapk pathways.

Gene networks and haloperidol-induced catalepsy.

The current study examined the changes in striatal gene network structure induced by short-term selective breeding from a heterogeneous stock for haloperidol response. Brain (striatum) gene expression data were obtained using the Illumina WG 8.2 array, and the datasets from responding and non-responding selected lines were independently interrogated using a weighted gene coexpression network analysis (WGCNA). We detected several gene modules (groups of coexpressed genes) in each dataset; the membership of the modules was found to be largely concordant, and a consensus network was constructed. Further validation of the network topology showed that using approximately 35 samples is sufficient to reliably infer the transcriptome network. An in-depth analysis showed significant changes in network structure and gene connectivity associated with the selected lines; these changes were validated using a bootstrapping procedure. The most dramatic changes were associated with a gene module richly annotated with neurobehavioral traits. The changes in network connectivity were concentrated in the links between this module and the rest of the network, in addition to changes within the module; this observation is consistent with recent results in protein and metabolic networks. These results suggest that a network-based strategy will help identify the genetic factors associated with haloperidol response.